arXiv:0708.1894 



Cosmological Constraints on New Agegraphic Dark Energy 

Hao WeiB 

Department of Physics and Tsinghua Center for Astrophysics, 
Tsinghua University, Beijing 100084, China 

Rong-Gen CaiQ 

Institute of Theoretical Physics, Chinese Academy of Sciences, P.O. Box 2735, Beijing 100080, China 



ABSTRACT 

In this work, we consider the cosmological constraints on the new agegraphic dark energy (NADE) 
proposed in arXiv:0708.0884 by using the observational data of type la supernovae (SNIa), the 
shift parameter from cosmic microwave background (CMB) and the baryon acoustic oscillation 
(BAO) peak from large scale structures (LSS). Thanks to its special analytic features in the 
radiation-dominated and matter-dominated epochs, NADE is a single-parameter model in prac- 
tice because once the single model parameter n is given, all other physical quantities of NADE 
can be determined correspondingly. The joint analysis gives the best-fit value (with la uncer- 
tainty) n = 2.716t£;iJg, and the derived fi m o, f2 g o and w q o (with la uncertainties) are 0.295^o 020> 
0.705^0 020 an d — 0-794lo'oo5 , respectively. In addition, we find that the coincidence problem could 
be solved naturally in the NADE model provided that n is of order unity. 
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I. INTRODUCTION 



In general relativity, one can measure the spacetime without any limit of accuracy. In quantum me- 
chanics, however, the well-known Heisenberg uncertainty relation puts a limit of accuracy in these mea- 
surements. Following the line of quantum fluctuations of spacetime, Karolyhazy and his collaborators [l[ 
(see also Q) made an interesting observation concerning the distance measurement for Minkowski space- 
time through a light-clock Gedanken experiment, namely, the distance t in Minkowski spacetime cannot 
be known to a better accuracy than 

st = \tyh i '\ (i) 

where A is a dimensionless constant of order unity 0, 0]. We use the units h — c = ks = 1 throughout 
this work. Thus, one can use the terms like length and time interchangeably, whereas l p — t p = l/m p 
with l p , t p and m p being the reduced Planck length, time and mass, respectively. 

The Karolyhazy relation |T]) together with the time-energy uncertainty relation enables one to estimate 
an energy density of the metric quantum fluctuations of Minkowski spacetime 0, H| . Following 0, Q , 
with respect to Eq. ([1]) a length scale t can be known with a maximum precision St determining thereby 
a minimal detectable cell St 3 ~ t^t over a spatial region t 3 . Such a cell represents a minimal detectable 
unit of spacetime over a given length scale t. If the age of the Minkowski spacetime is t, then over a 
spatial region with linear size t (determining the maximal observable patch) there exists a minimal cell 
St 3 the energy of which due to time-energy uncertainty relation can not be smaller than 0, Q 

S it »~t- 1 . (2) 

Therefore, the energy density of metric fluctuations of Minkowski spacetime is given by @, Q 

Est* 1 ^ 
Pq ~ St 3 ~ t*t* ~ ' 1 j 

We refer to the original papers [2, [1[ for more details. It is worth noting that in fact, the Karolyhazy 
relation (p} and the corresponding energy density ([3]) have been independently rediscovered later for 
many times in the literature (see e.g. [1, [1, Hi ) ■ 

Based on the energy density a so-called agegraphic dark energy (ADE) model was proposed in 
There, as the most natural choice, the time scale t in Eq. ((3]) is chosen to be the age of the universe 

where a is the scale factor of our universe; H = a/a is the Hubble parameter; an over dot denotes the 
derivative with respect to cosmic time. Then, the ADE model was extended to include the interaction 
between ADE and background matter The statefinder diagnostic and w — w' analysis for the ADE 
models without and with interaction were performed in [§]. The ADE was constrained by using some 
old high redshift objects [l(| and type la supernovae (SNIa), the shift parameter from cosmic microwave 
background (CMB) and the baryon acoustic oscillation (BAO) peak from large scale structures (LSS) [ill ]. 

On the other hand, a new model of ADE was proposed in [12], where the time scale in Eq. @ is chosen 
to be the conformal time r\ instead of the age of the universe. This new agegraphic dark energy (NADE) 
has some new features different from the ADE proposed in Q . We will briefly review the main points of 
NADE model in the next section. 

Thanks to its special analytic features in the radiation-dominated and matter-dominated epochs, NADE 
is a single-parameter model in practice, unlike the two-parameters ADE model [Til ] (this point will be 
explained below). If the single model parameter n is given, all other physical quantities of NADE are 
determined correspondingly. In this work, we find that the coincidence problem could be solved naturally 
in the NADE model provided that n is of order unity. In addition, we constrain the NADE by using 
the cosmological observations of SNIa, CMB and LSS. Different from the cosmological constraints on 
ADE [ll| , the single parameter n of NADE model can be strictly constrained. Then, the other physical 
quantities of NADE can also be determined tightly. 
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II. THE MODEL OF NEW AGEGRAPHIC DARK ENERGY 



A. Main points of the NADE model 

In the NADE model we choose the time scale in Eq. ^ to be the conformal time rj instead, which 
is defined by dt = adrj [where t is the cosmic time, do not confuse it with the t in Eq. JHJ)]. Therefore, 
the energy density of NADE reads 

3n 2 m 2 

where the numerical factor 3n 2 is introduced to parameterize some uncertainties, such as the species of 
quantum fields in the universe, the effect of curved spacetime [since the energy density in Eq. ([3]) is derived 
for Minkowski spacetime], and so on. Since both the numerical factor A in Eq. (|TJ) and the coefficient 
in time-energy uncertainty relation Eq. ^ are of order unity, it is anticipated that the parameter n in 
Eq. ([5|) which comes from Eqs. ([1]) and ([2]) is also of order unity. The conformal time i] is given by 

/dt f da 

If we write 77 to be a definite integral, there will be an integration constant in addition. Note that fj = l/a. 

In this work, we consider a flat Friedmann-Robertson- Walker (FRW) universe containing NADE and 
pressureless matter. Note that the assumption of flatness is motivated by the inflation scenario and is 
also consistent with the observation of CMB 26|. The corresponding Friedmann equation reads 

H2 = ^2 (Pm + Pi) > (?) 



where p m is the energy density of pressureless matter. It is convenient to introduce the fractional energy 
densities Qi = pi / '(3m 2 H 2 ) for i — m and q. From Eq. (|5|), the corresponding fractional energy density 
of NADE is given by 

By using Eqs. ((5]) — (JSj) and the energy conservation equation p m + 3Hp rn — 0, we find that the equation 
of motion for Q, q is given by [12] 

dtt q _ Q q n n s f o 2 ^/^^ q 



3 __v_« . (9) 
da a y n a I 

From the energy conservation equation p q + 3H(p q + p q ) = 0, as well as Eqs. ([8]) and ([5]), it is easy to 
find that the equation-of-state parameter (EoS) of NADE w q = p q / p q is given by [HJ 



Wq = -1 + ^-^. (10) 
3n a 

Obviously, the scale factor a enters Eqs. © and (fTTJ)) explicitly. When a — > 00, Q q — * 1, thus w q — > — 1 
in the late time. When a — * 0, Vl q — + 0, we cannot obtain w q from Eq. (fTUj) directly. Let us consider the 
matter-dominated epoch, H 2 oc p m oc a -3 . Thus, a 1 / 2 da oc dt = adrj. Therefore, r\ oc a 1 / 2 . From Eq. [[5]). 
p q oc a~ 1 . From the energy conservation equation p q + 3Hp q (l + w q ) = 0, we obtain that w q = —2/3 in 
the matter-dominated epoch. Since p m oc a~ 3 , it is expected that fl q oc a 2 . Comparing w q = —2/3 with 
Eq. (HUl), we find that tt q = n 2 a 2 /4 in the matter-dominated epoch as expected. For a <C 1, if n is of 
order unity, f2 g <C 1 naturally. On the other hand, one can check that £l q — n 2 a 2 /4 satisfies 

dn q = n q ( 2 y/n;\ 

rfa ay n a J ' 
which is the approximation of Eq. © for 1 — Q g ~ 1. Therefore, all things are consistent. 
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In fact, we can further extend our discussion to include the radiation-dominated epoch [121 ]. However, 
since our main aim of this work is to constrain NADE by using the cosmological observations of SNIa, 
CMB and LSS which only concern the matter-dominated epoch, here we skip the radiation-dominated 
epoch. Instead, we only give a brief summary of the results obtained in [l2l ]. namely, in the radiation- 
dominated epoch, w q = —1/3 whereas tt q — n 2 a 2 ; in the matter-dominated epoch, w q = —2/3 whereas 
ftq = n 2 a 2 /4; eventually, the NADE dominates; in the late time w q — > — 1 when a — > oo, the NADE 
mimics a cosmological constant. We refer to the original paper [l§] for more details. 




FIG. 1: The resulting Q. q {z), Q. m [z) = 1 — Q q (z) and w q (z) obtained from Eqs. pip and (|10[l with the initial 
condition f2 g (zim) = n 2 (l + Zi n i)~ 2 /4 at Zi n i = 2000 for different n. 



5 



B. NADE and the coincidence problem 

It is worth noting that both fl q = n 2 a 2 and n 2 a 2 /4 <C 1 naturally in radiation-dominated and matter- 
dominated epochs where a < 1. These results are obtained without any additional assumption. It is 
natural to consider the coincidence problem. In many dark energy models, especially the ACDM model, 
Qde should be fine-tuned to be extremely small in the early time, in order to interpret the fact that 
the energy densities of dark energy and matter are comparable today. In some dynamical dark energy 
models, the tracker behavior or attractors in the interacting dark energy models are used to alleviate the 
coincidence problem. 

Here, we see that the NADE model could shed new light on the coincidence problem in a different way. 
The key point is that as mentioned above tt q of NADE can be obtained analytically in both radiation- 
dominated and matter-dominated epochs. In particular, fl q = n 2 et 2 /4 = n 2 (l + z)~ 2 /4 in the matter- 
dominated epoch, where z = a -1 — 1 is the redshift (we set ao = 1 throughout; the subscript "0" indicates 
the present value of the corresponding quantity). Therefore, we can use ^(z^) = n 2 (l + Zi n i)~ 2 /4 at 
any Zi n i which is deep enough into the matter-dominated epoch to be the initial condition to solve the 
differential equation of il q . We rewrite Eq. @ as 

3(1 + Z)- 1 - -JSh . (11) 
n 

As well-known, the energy densities of matter and radiation is equal at redshift z eq which is given by 
1 + z eq ~ 2.32 x 10 4 n m0 h 2 pj|, where h is defined by the Hubble constant H = 100 hkm/s/Mpc. For 
instance, z eq ~ 3409 for Q m o = 0.3 and h = 0.7 whereas z eq ~ 2505 for O, m o — 0.3 and h — 0.6 [lj|. 
Therefore, we choose Zi n \ — 2000 which is deep enough into the matter-dominated epoch. In fact, one can 
check that if we use initial condition fi g (z, n ,) = n 2 (l + z.; n i)~ 2 /4 at Zi n i = 2000, the resulting fl q (z) from 
Eq. (TTTJ) for n = 3 relatively deviates from n 2 (l + z)~ 2 /4 less than approximately 1% in the large interval 
10 < z < 2000. Thus, if one uses other Zi n i, the results deviate from the ones of Zi n i = 2000 negligibly. 
This also justifies the choice of initial condition £l q (zi n i) = n 2 (l + Zi n i)~ 2 /4 at Zj n j = 2000. In Fig. [TJ 
we present the resulting fl q (z), fl m (z) = 1 — £l q (z) and w q (z) from Eqs. (fTT|) and (fT0|) for different n. 
We only plot the redshift range < z < 20 in Fig. [TJ since fl q ~ and w q ~ —2/3 closely in the range 
20 < z < 2000. From Fig. [TJ it is easy to see that £l q o and fl m Q = 1 — fl q o are comparable provided 
that n is of order unity. We stress that this result is obtained very naturally, without any additional 
assumption. Therefore, the coincidence problem could be solved naturally in the NADE model. 



_i = -n,(i-n,) 



III. CONSTRAINTS ON NADE FROM SNIA 



In this section, we constrain NADE by using the latest 182 SNIa Gold dataset [Ta|. The latest 182 
SNIa Gold dataset compiled in [l6[ provides the apparent magnitude m(z) of the supernovae at peak 
brightness after implementing corrections for galactic extinction, K-correction, and light curve width- 
luminosity correction. The resulting apparent magnitude m{z) is related to the luminosity distance 
cLl{z) through (see e.g. (l7l |) 

m th {z) = M(M, H ) + 5 log 10 D L {z), (12) 

where 

D L {z) = {l + z) f (13) 
Jo E(z;p) 

is the Hubble-free luminosity distance HoCLl/c in a spatially flat FRW universe (c is the speed of light); 
p denotes the model parameters; E(z) = H(z)/Ho; and 



/ c tj— i " 

M = M + 5 log 10 f ) + 25 = M - 5 log 10 h + 42.38 (14) 



G 



is the magnitude zero offset; h is Ho in units of 100 km/s/Mpc; the absolute magnitude M is assumed to 
be constant after the corrections mentioned above. The data points of the latest 182 SNIa Gold dataset 
compiled in [1 61 ] are given in terms of the distance modulus 

^obs(zi) = m obs (zi) - M. (15) 

On the other hand, the theoretical distance modulus is defined as 

m.h{z%) = m th (z. t ) - M = 51og 10 D L {zi) + fi , (16) 

where 

/z =42.38-51og 10 fc. (17) 
The theoretical model parameters are determined by minimizing 

2 / \ \ \P>obs(Zi) - IJ-thjzi)] 2 , . 

Xsn(P)=X, ^ > ( 18 ) 

where a is the corresponding la error; p denotes the model parameters. The parameter /io is a nuisance 
parameter but it is independent of the data points. One can perform an uniform marginalization over 
Hq. However, there is an alternative way. Following (l7l [Til. [l9| . the minimization with respect to /io can 
be made by expanding the Xsn oi Eq. (O with respect to /x as 

X 2 SN (p) =A-2noB + i$C ) (19) 

where 

is2 [m obs (zi) - mth(z t \ fi = 0, p)f 



^m obs (z i )~mth(z i - 1 fio = 0,p) 

(p) = ^ ^ — jv) ' 



c = 



182 

£ 

8 = 1 



.(*) 



Eq. (fT9|) has a minimum for /io = B/C at 

X 2 SJV (p)=A(p)-^. (20) 

Therefore, we can instead minimize Xsn which is independent of /io, since Xsn min = Xsn min obviously. 
It is worth noting that the corresponding h can be determined by /io — B/C for the best-fit parameters. 

In the ADE model [111 ], there are two independent model parameters, namely, n and il m o- To obtain 
flq(z) which is critical in E(z), one should solve the corresponding differential equation of £l q with the 
initial condition Q q (z — 0) = fl q o = 1 — fimo, similar to the case of holographic dark energy model (see 
e.g. [23|). In the NADE model, however, the situation is different. As mentioned in Sec. Ill Bl if n is given, 
we can obtain tt q (z) from Eq. (fTTj) with the initial condition Q q (zi n i) = n 2 (l + Zi n i)~ 2 /4 at Zi n i = 2000 (or 
any Zi n i which is deep enough into the matter-dominated epoch), instead of Cl q (z = 0) = Cl q o = 1 — fl m o- 
Then, all other physical quantities, such as fl m (z) = 1 — Q, q (z) and w q (z) in Eq. (fTO]) . can be obtained 
correspondingly. So, £l m o = ^m(z = 0) and Cl q o — Q q (z — 0) are not independent model parameters. The 
only model parameter is n. Therefore, the NADE model is a single-parameter model in practice, unlike 
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the two-parameters ADE model [llj. To our knowledge, it is the third single-parameter cosmological 
model besides the well-known ACDM model and the DGP braneworld model [33| . From Eq. ([7]), we have 



E[z) = 



n m0 (i + zf 



-,1/2 



i - n q (z) 



(21) 



If the single model parameter n is given, we can obtain Q. q (z) from Eq. pip . And then, we get f2. 



mO 



1 — = 0). Therefore, E(z) is in hand. So, we can find the corresponding Xsn- In Fig. El we present 
the x 2 — Xsn an d the corresponding likelihood C oc e~ x / 2 for < n < 10. We find that the best-fit 
model parameter is n = 2.954 while Xmin = 160.255. The corresponding h — 0.623. In Table [H we 
present the best- fit value of n and the derived f2 m o, Qqo and w q o with la and 2a uncertainties. On the 
other hand, we also fit the ACDM model to the same SNIa dataset and find that Xmin, A = 158.749 for 
the best- fit parameter S!^ = 0.342 (the corresponding h = 0.626). So, fitting to the latest 182 SNIa 
Gold dataset [II], the ACDM model is slightly better than the NADE model. 




FIG. 2: The x 2 = Xsn an d the corresponding likelihood for < n < 10. These results are obtained from the 
latest 182 SNIa Gold dataset compiled in [la ]. 



Uncertainty 


n 






W q0 


la 
2a 


2 Q54+0-264 
z ' MO4 -0.245 

2.954tg;g» 


u.zoo_ 038 

n oc:c+0.086 
u.zoo_ 072 


n 74^+0-038 

u. i4o_ 041 

n 74c;+0 072 
u - '*°-0.086 


-0.805l°;°l 2 2 

-0.80518-™ 



TABLE I: The best-fit value of n and the derived Q m o, QqO and w a o with lcr and 2a uncertainties. These results 
are obtained from the latest 182 SNIa Gold dataset compiled in [la] . 
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n 



FIG. 3: The \ 2 an( A the corresponding likelihood for < n < 10. These results are obtained from the combined 
SNIa, CMB and LSS data. 



IV. CONSTRAINTS ON NADE FROM SNIA, CMB AND LSS 



Although the constraints on NADE from SNIa alone are fairly tight as shown in Table HI we can 
further make the constraints tighter by combining SNIa with other complementary observations. In the 
literature, the shift parameter R from CMB [20l. T2ll] and the parameter A of BAO measurement from 
LSS [H[ are used extensively, see [HlQilHII for examples. It is commonly believed that both R and A are 
model-independent and contain the essential information of the full CMB and LSS BAO data (however, 
see also e.g. [24( and |25|). Notice that both R and A are independent of Hubble constant H . In fact, 
the shift parameter R is defined by [13, HH 

where z rec = 1089 is the redshift of recombination. The shift parameter R relates the angular diameter 
distance to the last scattering surface, the comoving size of the sound horizon at z re c and the angular 
scale of the first acoustic peak in CMB power spectrum of temperature fluctuations [13, Hl| . The value of 
R is determined to be 1.70 ± 0.03 [Hj] from the WMAP 3-year (WMAP3) data On the other hand, 
the parameter A of the measurement of BAO peak in the distribution of SDSS luminous red galaxies [27| 
is given by [H| 

A ^ n^E(z„r 1/3 

where z b = 0.35. In 22], the value of A is determined to be 0.469 (n s /0.98)~ - 35 ± 0.017, here the scalar 
spectral index n s is taken to be 0.95 from the WMAP3 data [261 ] - 

We perform a joint analysis of SNIa, the shift parameter R from CMB and the BAO peak measurement 
A from LSS data to constrain the NADE model. The total \ 2 is given by 

X 2 = Xsn + Xcmb + xtss> ( 24 ) 

where X 2 SN is given in Eq. O, Xcmb = ( R - Robs) 2 /v\ and xlss = ( A ~ A obs ) 2 / (J 2 A . In Fig. [3 we 
present the \ 2 an( i the corresponding likelihood C oc e~ x / 2 for < n < 10. We find that the best-fit 



1 f Zb dz 
Zb Jo E(z) 



(23) 
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Uncertainty 


n 




n g o 


WqO 


la 
2a 


2 71 6+ - 111 

' ±u -0.109 


u.zyo_ 020 
295+ ' 041 

U ' ZJJ -I) .038 


U. /UO_ 020 
U - ' UO -0.041 


U - ' y4 -0.005 

-0 704+0 11 

U- '^-0.010 



TABLE II: The best-fit value of n and the derived flmo, QqO and w q o with l<r and 2a uncertainties. These results 
are obtained from the combined SNIa, CMB and LSS data. 

model parameter is n = 2.716 while x m in = 161.467. The corresponding h = 0.616. We present the 
best-fit value of n and the derived fl m o, QqO and w q o with la and 2a uncertainties in Table [TT1 Obviously, 
these constraints in Table HI1 are strict enough. On the other hand, we also fit the ACDM model to the 
same combined SNIa, CMB and LSS data and find that Xmin A = 162.886 for the best-fit parameter 
fl^ = 0.288 (the corresponding h = 0.637). So, fitting to the combined SNIa, CMB and LSS data, the 
NADE model is slightly better than the ACDM model. This makes the NADE model more attractive. 



V. CONCLUSION 

In this work, we consider the NADE model proposed in fl2j |. which might have some physical motiva- 
tions connected to the quantum fluctuations of spacetime. Thanks to its special analytic features in the 
radiation-dominated and matter-dominated epochs, NADE is a single-parameter model in practice, unlike 
the two-parameters ADE model 11]. If the single model parameter n is given, all other physical quan- 
tities of NADE can be determined correspondingly. To our knowledge, it is the third sin gle-p arameter 
cosmological model besides the well-known ACDM model and the DGP braneworld model [33 |. We find 
that the coincidence problem could be solved naturally in the NADE model provided that the single 
model parameter n is of order unity. In addition, we constrain NADE by using the cosmological obser- 
vations of SNIa, CMB and LSS. The joint analysis gives the best-fit parameter (with la uncertainty) 
n = 2.716±g iog. The derived ft m0 , tt q0 and w q0 (with la uncertainties) are 0.2951^20 , 0.705;^;^o 
and —0.794^0 005! respectively. These constraints are strict enough. On the other hand, fitting to the 
combined SNIa, CMB and LSS data, we find that the NADE model is slightly better than the ACDM 
model. This makes the NADE model more attractive. 

If we add other observations such as Chandra X-ray observations 28], observational H(z) data [29| . 
the lookback time data compiled in [3(| and so on, it is expected that the constraints on NADE model 
will be very tight and the resulting parameters can be used to make some exact predictions of NADE. 
Of course, the other SNIa datasets such as SNLS [3l[ and ESSENCE [33] are also useful to constrain the 
NADE model. It is worth noting that in this work we only used the parameters R and A from WMAP3 
and SDSS data, the resulting constraint on n is primary in some sense, and therefore it is of great interest 
to constrain the NADE model by using the global fitting to the full CMB and LSS data via Markov Chain 
Monte Carlo (MCMC) analysis. This is an issue which deserves further investigation in the future. 



ACKNOWLEDGMENTS 



We thank the anonymous referees for quite useful comments and suggestions, which help us to improve 
this work. We are grateful to Prof. Shuang Nan Zhang for helpful discussions. We also thank Minzi Feng, 
as well as Hui Li, Yi Zhang, Xing Wu, Li-Ming Cao, Xin Zhang, Jian Wang and Bin Hu, for kind 
help and discussions. This work was supported in part by a grant from China Postdoctoral Science 
Foundation, a grant from Chinese Academy of Sciences (No. KJCX3-SYW-N2), and by NSFC under 
grants No. 10325525, No. 10525060 and No. 90403029. 



10 



[1] F. Karolyhazy, Nuovo Cim. A 42, 390 (1966); 

F. Karolyhazy, A. Frenkel and B. Lukacs, in Physics as Natural Philosophy, edited by A. Simony and 
H. Feschbach, MIT Press, Cambridge, MA (1982); 

F. Karolyhazy, A. Frenkel and B. Lukacs, in Quantum Concepts in Space and Time, edited by R. Penrose 

and C. J. Isham, Clarendon Press, Oxford (1986). 
[2] M. Maziashvili, Int. J. Mod. Phys. D 16, 1531 (2007) [gr-qc/0612lT(j] . 
[3] M. Maziashvili, Phys. Lett. B 652, 165 (2007) [arXiv:070 5.0924 . 
[4] N. Sasakura, Prog. Theor. Phys. 102, 169 (1999) [hep-th/9903146| . 
[5] Y. J. Ng and H. Van Dam, Mod. Phys. Lett. A 9, 335 (1994); 

Y. J. Ng and H. Van Dam, Mod. Phys. Lett. A 10, 2801 (1995); 

S. Lloyd and Y. J. Ng, Sci. Am. 291, 52 (2004). 



W. A. Christiansen, Y. J. Ng and H. van Dam, Phys. Rev. Lett. 96, 051301 (2006) gr-qc/0508121 
M. Arzano, T. W. Kephart and Y. J. Ng, Phys. Lett. B 649, 243 (2007) [gr-qc/0605lT7] ; 



Y. J. Ng, Phys. Lett. B 657, 10 (2007) gr-qc/0703096 



[7] R. G. Cai, Phys. Lett. B 657, 228 (2007) |arXiv:0707. 4049 . 

[8] H. Wei and R. G. Cai, larXiv:0707.4052l [hep-th] . 

[9] H. Wei and R. G. Cai, Phys. Lett. B 655, 1 (2007) [arXiv:0707.4526] . 
[10] Y. Zhang, H. Li, X. Wu, H. Wei and R. G. Cai. larXiv:0708.1214l [astro-ph], 
[11] X. Wu, Y. Zhang, H. Li, R. G. Cai and Z. H. Zhu, larXiv:0708.0349 l [astro-ph]. 
[12] H. Wei and R. G. Cai, Phys. Lett. B 660, 113 (2008) [arXiv:0708.0884] . 

[13] E. W. Kolb and M. S. Turner, The Early Universe, Addison- Wesley Publishing Company, USA (1990) [in 

the series Front. Phys. 69, 1 (1990)]. 
[14] W. L. Freedman et al, Astrophys. J. 553, 47 (2001) [astro-ph/0012376l . 
[15] A. Sandage, G. A. Tammann, A. Saha, B. Reindl, F. D. Macchetto and N. Panagia, 

Astrophys. J. 653, 843 (2006) [astro- ph/0603647] . 
[16] A. G. Riess et al. [Supernova Search Team Collaboration], astro-ph/0611572 

The numerical data of the full sample are available at 

http://braeburn.pha.jhu.edu/~ariess/R06 or upon request to ariess@stsci.edu 



[17] S. Nesseris and L. Perivolaropoulos, Phys. Rev. D 72, 123519 (2005) astro-ph/0511040 
[18] L. Perivolaropoulos, Phys. Rev. D 71, 063503 (2005) [astro-ph/0412308l. 



[19] E. Di Pietro and J. F. Claeskens, Mon. Not. Roy. Astron. Soc. 341, 1299 (2003) astro-ph/0207332 



[20] J. R. Bond, G. Efstathiou and M. Tegmark, Mon. Not. Roy. Astron. Soc. 291, L33 (1997) astro-ph/9702100 . 

[21] Y. Wang and P. Mukherjee, Astrophys. J. 650, 1 (2006) [astro-ph/06040511 . 

[22] D. J. Eisenstein et al. [SDSS Collaboration], Astrophys. J. 633, 560 (2005) [astro-ph/0501171| . 

[23] X. Zhang and F. Q. Wu, Phys. Rev. D 76, 023502 (2007) [astro-ph/0701405] . 

[24] O. Elgaroy and T. Multamaki, astro-ph/0702343| 

Y. Wang and P. Mukherjee, Phys. Rev. D 76, 103533 (2007) [astro-ph/0703780| . 

[25] S. Rydbeck, M. Fairbairn and A. Goobar, JCAP 0705, 003 (2007) [astro-ph/0701495| . 



[26] D. N. Spergel et al. [WMAP Collaboration], Astrophys. J. Suppl. 170, 377 (2007) astro- ph/0603449 



L. Page et al. [WMAP Collaboration], Astrophys. J. Suppl. 170, 335 (2007) astro-ph/0603450 ; 
G. Hinshaw et al. [WMAP Collaboration], Astrophys. J. Suppl. 170, 288 (2007) [astro-ph/06 03451 
N. Jarosik et al. [WMAP Collaboration], Astrophys. J. Suppl. 170, 263 (2007) astro-ph/0603452; . 
[27] M. Tegmark et al. [SDSS Collaboration], Phys. Rev. D 69, 103501 (2004) [astro-ph/0310723]; 



M. Tegmark et al. [SDSS Collaboration], Astrophys. J. 606, 702 (2004) astro-ph/0310725 



U. Seljak et al, Phys. Rev. D 71, 103515 (2005) astro-ph/0407372 ; 

J. K. Adelman-McCarthy et al. [SDSS Collaboration], Astrophys. J. Suppl. 162, 38 (2006) [astro-ph/0507711| ; 
K. Abazajian et al. [SDSS Collaboration], |astro-ph/0410239| |astro-ph/0403325| |astro-ph/0305492| 
M. Tegmark et al. [SDSS Collaboration], Phys. Rev. D 74, 123507 (2006) astro-ph/0608632]. 



11 



[28] S. W. Allen et al., Mon. Not. Roy. Astron. Soc. 353, 457 (2004) [astro-ph/04"0 5340 ; 

S. W. Allen et al. . llarXiv:0706.0033 1 [astro-ph] . 
[29] R. Jimenez, L. Verde, T. Treu and D. Stern, Astrophys. J. 593, 622 (2003) astro- ph/0302560 ; 

J. Simon, L. Verde and R. Jimenez, Phys. Rev. D 71, 123001 (2005) [astro-ph/0412269l ; 

L. Samushia and B. Ratra, Astrophys. J. 650, L5 (2006) [astro-ph/0607301| ; 

H. Wei and S. N. Zhang, Phys. Lett. B 644, 7 (2007) astro-ph/0609597 ; 

H. Wei and S. N. Zhang, Phys. Lett. B 654, 139 (2007) [arXiv:0704.3330] . 
[30] S. Capozziello, V. F. Cardone, M. Funaro and S. Andreon, Phys. Rev. D 70, 123501 (2004) [astro-ph/0410268| . 
[31] P. Astier et al. [SNLS Collaboration], Astron. Astrophys. 447, 31 (2006) [astro-ph/0510447| . 
[32] W. M. Wood-Vasey et al. [ESSENCE Collaboration], Astrophys. J. 666, 694 (2007) [astro-ph/0701041] ; 

G. Miknaitis et al. [ESSENCE Collaboration], Astrophys. J. 666, 674 (2007) [astro-ph/0701043] . 
[33] G. R. Dvali, G. Gabadadze and M. Porrati, Phys. Lett. B 485, 208 (2000) hep-th/0005016 ; 

C. Deffayet, G. R. Dvali and G. Gabadadze, Phys. Rev. D 65, 044023 (2002) [astro-ph/0i 05068 . 



